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Instability of point defects 
in a two-dimensional nematic liquid crystal model 

Radii Ignat*, Luc Nguyenf Valeriy Slastikov* and Arghir Zarnesci# 


Abstract 

We study a class of symmetric critical points in a variational 2D Landau - de Gennes 
model where the state of nematic liquid crystals is described by symmetric traceless 
3x3 matrices. These critical points play the role of topological point defects carrying a 
degree | for a nonzero integer k. We prove existence and study the qualitative behavior 
of these symmetric solutions. Our main result is the instability of critical points when 
k ±1,0. 

1 Introduction 

1.1 Physical motivation 

The defining feature of nematic liquid crystals is the local orientational ordering of the 
molecules. Its main macroscopic manifestation is the emergence of certain patterns, called 
defects (points, lines or surfaces) where the local ordering, either disappears or changes 
abruptly. Defects determine a number of the most important features of liquid crystals, 
underlying spectacular phenomena and new prospective technologies, e.g. knotted disincli¬ 
nation lines, bistable displays, control of nanoparticle suspensions (see my These defects 
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are often analysed in comparison with topological singular phenomena appearing in other 
fields of condensed and soft matter physics, such as superconductivity, materials science, 
physics of polymers and even cosmology. 

There exist several competing continuum liquid crystal theories describing the local orien¬ 
tational ordering by a specific order parameter (see The most comprehensive 

and widely accepted continuum theory of nematic liquid crystals is the Landau-de Gennes 
theory [8]. It uses as an order parameter the so-called Q-tensor (a traceless, symmetric 3x3 
matrix) so that the analysis is carried out in the five-dimensional space J^ 0 °f Q-tensors: 


y»= Qe K 3X3 : Q = Q t i tr(Q) = 0 


= < s ( n <g> n - I 3 ) + r ( m <g> m - I 3 I : s, r £ M., n, m £ S 2 , n ■ m = 0 >, 


where § 2 is the unit sphere in M 3 , I 3 is the 3x3 identity matrix and (n <g) nj.. = riiTij for 
1 < i,j < 3. 

The central object in the Landau-de Gennes theory is the free energy functional <^(Q); 
in fact, stable equilibrium configurations of the liquid crystalline system in G C (d = 2, 3) 
correspond to local minimisers of Landau-de Gennes energy. In the simplest form, the free 
energy of a liquid crystal is given by 


&(Q) 


|VQ| 2 + f bulk (Q) dx, Q G Hl oc (n,Y 0 ). 


( 1 . 1 ) 


The simplest bulk potential fbuik(Q ) that captures the main physical characteristics is taken 
to be of the form 

fbuik(Q) = -ytr (Q 2 ) - jtr(Q 3 ) + j (tr (Q 2 )f , 

where a 2 , b 2 , c 2 > 0 are material constants. Note that the minimum set of the bulk potential 
fbuik(Q) is given by the set of uniaxial Q tensors (i.e., Q has two equal eigenvalues): 


= js+ yn®n- -I 3 

with the constant order parameter s + given by 


n 6 


b 2 + v& 4 + 24a 2 c 2 „ 

S+ = - 4 ?- >0 ' 


( 1 . 2 ) 


(1.3) 


The critical points of the energy functional bP(Q) satisfy the Euler-Lagrange equation: 


A Q = —a 2 Q - b 2 [Q 2 --\Q\ 2 I 3 } + c 2 Q\Q\ 2 in G, 

O 


(1.4) 
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where | \Q\ 2 = |tr(Q 2 ) is the Lagrange multiplier associated to the traceless constraint. It is 
known that any 77^(12, J^o)-solution of (11.4(1 is smooth in 12. The solutions of (II .4(1 describe 
the defects patterns, the simplest and most common being the point defects (see [6], 22, |23]). 
The analytical investigation of their structure and profile generates very challenging nonlinear 
analysis problems. 

The goal of this article is to investigate the profile and stability properties of point defects 
appearing for a certain type of symmetric solutions of (11.41) in the two dimensional case 

12 = M 2 . 


The boundary conditions imposed for these solutions are taken to be: 


Q(x) —y Qk(x) := s + ( n(x) <g) n(x) — -/ 3 ) as |x| —> oo, 


where the map n : 12 —> S 2 is given in the polar coordinates by 


(1.5) 


n(x) — (cos(|</?), sin(|</p), 0) , r > 0, (p G [0, 27t), (1.6) 

where k 6 Z and x = (rcostp, rsintp). Note that if we consider Qas an MP 1 -valued map 
on IR 2 \ {0}, then it has degree kj 2 about the origin. (For a definition of the degree for MP 1 - 
valued maps, see for instance (4j, p.685 — 686). This model can be seen as the 2D reduction 
of the physical situation of a 3D cylindrical boundary domain, with so-called “homeotropic” 
boundary conditions where the configurations are invariant in the vertical direction (see for 
instance 0 )- 


1.2 The /c-radially symmetric solutions 

We will focus on the following type of symmetric solutions of (11.41) in the two-dimensional 
domain 12 = M 2 that carry a topological information through the boundary condition (11.51) . 

Definition 1.1. For k e Z \ {0}, we say that a Lehesgue measurable map Q : 12 — > J^o is 
k-radially symmetric if the following conditions hold for almost every x = (ay, x 2 ) G 12: 


(HI) The vector e 3 = (0,0,1) is an eigenvector of Q(x). 

(H2) The following identity holds 

Q ^P 2 (Pfc(V')5)^ = 7 Zk{if)Q(x)TZ t k {'if), for almost every if el, 

where x = (ay, £2,0), P 2 : K 3 —> M 2 is the projection given as P 2 (ay, x 2 , x 3 ) = (xi,x 2 ) and 


Tlkipf) 


cos (| if) 

sin(fV’) 

0 


— sin(|-0) 

cos(|^) 
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0 

0 

1 


(1.7) 









is the winding rotation around the vertical axis e 3 . 

Remark 1.2. If k is an odd integer, then a map Q G satisfying (H2) automati¬ 

cally verifies (HI) (see Proposition 1 2. II) . 


We will show that the fc-radially symmetric solutions of (11.411 have a simple structure: 

Proposition 1.3. Let \ {0}. If Q G ^ 0 ) a k-radially symmetric solution 

of the Euler-Lagrange equations (jl .4[) on Q = M. 2 satisfying the boundary conditions (11.51) , 
then Q is smooth and has the following form for every 1 6 R 2 ; 


Q(x) = u( \x\)V2 ^n(ar) ®n(x) - +v(\x\)^ ^e 3 ® e 3 - ^/ 3 ^ , (1.8) 

u>/iere n is given in (jl.6jl . / 2 = / 3 — e 3 <g) e 3 , u G C 2 ([0, 00)) D C°°((0, 00 )), v G C°°([0, 00)) 
and u and v satisfy the following system of ODEs in (0, 00): 


\v" + ^ ^ =U 

— a 2 + y" (b 2 v + c 2 (u 2 + v 2 ) 

[v" + y =v 

—a 2 — b 2 v + c 2 (u 2 + v 2 ) 


V6 


subject to boundary conditions: 


M (0) = 0, v'(0) = 0, u(+ 00 ) = -J=s+, u(+oo) = --^=s+. 


(1.9) 


( 1 . 10 ) 


Conversely, if u G Hf oc ([0, 00); rdr) O L 2 OC ([0, 00); ^) and v G H) oc ([0, 00); rdr) satisfy (jl,9[) 
with the boundary condition u(+ 00) = -j^s + and v(+oo) = —^s + , then the tensor Q defined 
by (ELBD belongs to H) oc (M. 2 , J/q) and is a k-radially symmetric smooth solution of (11.41) - (1 1.5 f) . 


Analysing the above ODE system, we construct solutions of (1 1.9 [) - (II. 10p using variational 
methods that lead to fc-radially symmetric solutions of the Euler-Lagrange equations (11.41) 
with the boundary conditions fll.5|) . 


Theorem 1.4. Let a 2 , b 2 , c 2 > 0 be any fixed constants and k G Z\{0}. There exist k-radially 
symmetric solutions Q G /// oc (R 2 , .5^) of (jl.4[) - (II,5jl having the form (II.80 . Moreover, the 
corresponding profiles (■ u , v) in (II .8[) satisfy the ODE system (jl.9[) - (II . 10)1 together with 


u > 0 and v < 0 in (0, 00 ). 


Remark 1.5. The case b 2 = 0 was studied in m on bounded domains. They showed that on 
bounded domains, the ODE system has a unique solution under the assumption that u > 0 
and v < 0. However, for infinite domain, the condition b 2 > 0 is essential in Theorem \1. f : 
there is no solution to the ODE system on (0, 00) with b 2 = 0 which satisfies u > 0 and 
v < 0. See Appendix\A 1 


Open problem 1.6. Forb 2 > 0, does the ODE system fll.9D -f ll.10D have a unique solution? 
(See Proposition \3.5\ for a statement to this effect in a special case.) 
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1.3 Instability of /c-radially symmetric solutions for k (E Z \{0,±1} 

Our main result concerns the local instability for all fc-radially symmetric critical points of 
& subject to (11.51) when fceZ\{ 0, ±1}: 

Theorem 1.7. Let a 2 , b 2 , c 2 > 0 be any fixed constants and fceZ\{ 0, ±1}. Any k-radially 
symmetric critical point Q of (11.11) with 0 = M 2 satisfying the boundary conditions (11.5ft is 
locally unstable, i.e. there is a perturbation P e C £°(supported in a bounded disk 
Br, such that the second variation Jz?[Q](P) < 0, where 

&[Q]{P) = L {i |v(Q + iP)|2 + / “‘ w + tp) ~ 5 |ve|2 “ dx 

= J {l|Vi=| 2 - ^|P| 2 - b 2 tr(P 2 Q) + 7 (|Q| 2 |P| 2 + 2|tr(QP)| 2 ) } <fc. (1.11) 

Open problem 1.8. Is it true that k-radially symmetric solutions of (11.4ft in R 2 subject to 
(11.5ft are stable for k — ± 1? 

Remark 1.9. This instability behaviour is drastically different from the case b 2 = 0 on a 
bounded disk Br centered at the origin. In J2J/, it was shown that the functional & with 
a boundary condition similar to (11 . 5[) has a unique global minimizer in H L (BR,y 0 ), and 
furthermore that minimizer is k-radially symmetric. The deeper reason for this seems to be 
related to the different structure of the minimum set of the potential fbuik, which for b 2 = 0 
is a AD sphere while for b 2 0 is the 2D real projective plane. 

There have been numerous numerical and analytical studies of two-dimensional point 
defects in the Landau - de Gemies framework [U [5j [7, E, T3[ [15 : , T71 [2T, 25] (also in micro- 
magnetics see e.g. mm )• Let us briefly mention a few papers that are directly related 
to this work. Our motivation came from the recent paper [9] where global minimisers of 
Landau-de Gennes energy are investigated on finite two-dimensional balls in the extreme 
low-temperature regime ( b 2 = 0) under fc-radially symmetric homeothropic boundary condi¬ 
tions. The authors show that there exists a unique global minimizer of the Landau-de Gennes 
energy which is fc-radially symmetric and provide the description of the ground state profile 
of a point defect of index k/2 in terms of the system of two ordinary differential equations 
(see (HD). More general domains and boundary conditions were treated analytically (see 
mum) and numerically (see mi). In |1] the Landau-de Gennes energy was investigated 
in a restricted three dimensional space of Q-tensors. The authors showed that in the case 
of small elastic constant the minimizers of Landau-de Gennes energy exhibit behavior sim¬ 
ilar to those of Ginzburg-Landau energy [3j, namely for boundary conditions of degree k/2 
there are exactly k vortices of degree ±1/2. In [5j [15] the minimizers of the full Landau-de 
Gennes energy were studied under non-orientable boundary conditions (which in our setting 
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amounts to k being odd). It was shown that in the low temperature regime and in the case 
of small elastic constant the minimizer has only one vortex. 

The paper is organised as follows: in the next section we provide the basic properties 
of the fc-radially symmetric maps that we study on balls Br of radius R G (0, oo]. In 
Section 3 we investigate the ODE system (11.91) on bounded domains and prove certain fine 
qualitative properties of solutions that will be used later. In Section 4 we show the existence 
of a fc-radially symmetric solution on the whole M 2 and investigate its behaviour at infinity. 
Finally, in Section 5 we investigate the stability of /c-radially symmetric solutions and show 
Theorem 11.71 Several open questions are also stated, some of them will be addressed in a 
forthcoming paper. 


2 Basic aspects of /c-radially symmetric maps, k ^ 0 


In order to classify fc-radially symmetric maps on balls Br centered at the origin with 
R G (0, oo] and k 0 (see Definition 11.11 for D = B R ), we introduce some notation. We 
define {e;} 3 =1 to be the standard basis in M 3 and denote, for p G [0, 2n), 

n = n(p) = (cos(|</?),sin(|</?),0) , m = m{p) = (— sin(|</?), cos(|</?), 0) . 

We endow the space of Q-tensors with the scalar product 

Q ■ Q = tr (QQ) 

and for any p G [0, 2i r), we define the following orthonormal basis in 


E o — \ h; \ e 3 ® e 3 


E i = Ei ((p) = \[2 


\ l 


-h 

2 


n($n- ~ii I , E 2 = E 2 (<p) = -j= (■n O 
E 3 = —^=(ei <8) e 3 + e 3 <g) ei), E 4 = -j= (e 2 <S> e 3 + e 3 <g) e 2 ). 


m + m®n) 


Obviously, only E 3 and E 2 depend on tp and we have 


dEi 


= kE - 2 and 


dE 2 


= -kEi. 


( 2 . 1 ) 


( 2 . 2 ) 


(2,3) 


d<p dp 

We prove the following characterization of property (H2) for H 1 (Br,S > o)-maps. 

Proposition 2.1. Let R G (0, cxo), k ^ 0 and Q g Ef 1 (B R , 5? 0 ) he a map that satisfies (H2) 
in B R .Then Q can he represented for a.e. x = r(cosp,smp) G B R : 0 


1 In these decompositions of Q(x), the angle ip defining Ei and E 2 is given by the phase of x G Br. 
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1. If k is odd, then 


2 

Q = 

i =0 

where w 0 G H l {(0, R)-,rdr) and uui,w 2 G R)-,r dr) n L 2 ((0,i?); ^dr). 

2. If k is even, then 

k k k k 

Q = 2_^ w i( r )Ei + (u’( r ) cos -(p + w(r) sin -<p)E 3 + (—w(r) cos -92 + w(r) sin ~(p)E A , 

i =0 


uAere w 0 G i7 1 ((0, i?);r dr) and w,w,wi,w 2 G i/ 1 ((0, i?); r dr) H L 2 ((0,/?); ^ dr). 

Proof. Any Q G H 1 (Br, S^q) can be represented as 

4 

Q(x) = YMx)Ej, x G Br, 

i= 0 

with Wi — Q ■ Ei for i — 0, ..., 4 and 

k?i 2 = X >, 2 

i =0 

, o lT _ ,9 k 2 . n 9 . 2k, dw 2 dvj ], , . 

and |VQ| =5^|Vwi| +^(w 1 +w 2 ) + ^(w 1 —-w 2 —), (2.4) 

where we used (12.3|) . Now we compute for ip G R: 

Rk(ip)e 1 = n(ip),IZ k (ip)e 2 = m(ip),TZ k (ip)e 3 = e 3 , 

7l k (ip)n(<p) = n{p + ip),7l k {ip)m(<p) = m{p + ip),ip G [0, 27t), 

so that we have for a.e. £ = r(cos</?, sin</?) G 5/?: 

Rk{fp)Q(x) Rl(ip) = w 0 {x)E 0 + wi(x)Ei(p + ip) + w 2 (x)E 2 (<p + ip) 

k k k k 

+ 1^3(0;) (cos —ipE 3 + sin —ipE/f) + in 4 (a;)(— sin —ipE 3 + cos —ipE/f). 

Therefore, hypothesis (H2) is equivalent (in polar coordinates) with: 

Wi(r, <p + ip)= Wi(r , 92), i = 0,1, 2 

k k 

w 3 (r, <p + ip) = w 3 (r , <p) cos - 104(r, (p) sin -ip, 

fc k 

w 4 pr, ip + ip) = w 3 (r , <p) sin -ip + w 4 (r, <p) cos --0, 
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for a.e. r G (0 ,R),<p G (0,27r),0 G M. Therefore, we deduce that u>; are independent 
of the angular variable p for i = 0,1,2. Since Q G H L (B R , 5?q), we obtain that w 0 G 
i/ 1 ((0, £?); r dr), Wi G /^((O, £?); r dr) D L 2 ((0, R); ^ dr ), i = 1, 2 and Ws,w 4 G ^(Br). It 
remains to characterize and w 4 . Let r G (0, R) so that w 3 and w 4 are continuous on dB r . 
(This is true because w 3 ,w 4 G H 1 (dB r ) C C 0, ^(dB r ) for a.e. r G (0,i?).) Then the above 
equalities for w 3 and w 4 hold for every p G [0, 27t) and ^ G t. Setting w(r) = w^(r, 0) and 
w(r) = —w 4 (r, 0), we get that 

k k 

Wz{r, 0) = w(r) cos —0 + w(r) sin —0, 
k k 

w 4 (r, 0) = w{r) sin —0 — w(r) cos —0, 

for every 0 G R. If A; is odd, the continuity of the 27r-periodic functions W 3 (r, •) and w 4 {r, ■) 
implies that w(r) = w(r) = 0 for a.e. r G (0, R). If k is even, then 


w 2 {x) + w\(x) = w 2 (r) + w 2 {r) 


k 2 


| Vw 3 | 2 (a;) + \X7w 4 \ 2 (x) = ( w') 2 (r) + ( w') 2 (r) + — ~(w 2 (r) +w 2 (r)) 

4 r z 


for a.e. x = r(cosp, sin<^) G B R . The proof is now completed. □ 

As a consequence, we deduce the following characterization of A;-radiaIly symmetric maps 
defined on balls B R \ 

Corollary 2.2. Let R G (0, 00 ) and k 0 0. A map Q G H l (B R , ,R 0 ) is k-radially symmetric 
if and only if Q can be represented as 


2 

Q = w i( r ) E ^ x = r ( cos sin P) G B Ri (2.5) 

i=0 


where w 0 G H l ({f), R);rdr) and Wi G ££ 1 ((0, R); r dr) D L 2 ((0, i?); £ dr), i = 1, 2. Moreover, 
we have 


&(Q) 

2?r 




3(w 2 + w 2 )] rdr. 



rdr 


Proof. Assume that Q is fc-radially symmetric. If k is odd, (12.5ft follows directly from Propo¬ 
sition [2T] If k is even, by (HI), e 3 is an eigenvector of Q and so the functions w and w 
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obtained in Proposition 12.11 are zero almost everywhere in (0,-R). In either case, we have 
proved (12.51) . The converse implication is obvious. Now, we compute for Q given by (12. 5 1) : 


tr (Q 3 ) = K - 3 K 2 + w l)\- 

The expression of & immediately follows. □ 

We now provide the proof of Proposition 11.31 with the characterization of fc-radially sym¬ 
metric solutions of (ll.4p - (ll.5p . In fact, we will prove the result for arbitrary balls B R with 
R G (0, 00 ]. The existence of such solutions is postponed to the next two sections. 

Proposition 2.3. Let k 7 ^ 0 and R G (0, 00 ]. If Q G Hf oc (B R ,^ 0 ) be a k-radially symmetric 
solution of the Euler-Lagrange equations (11.41) for = B R that satisfies the homeotropic 
boundary condition 


Q(x) 


s + ( n(x) <g) n(x) - -J 3 


on dB R 


( 2 . 6 ) 


(with the convention (II.5p if R = 00 ), then Q is smooth and 


Q(x) = v(r)E 0 + u(r)Ei(<p) for every x = r(cosip, sin</?) G B R , 

u G C 2 ([0,f?)) n C°°((0,R)), v G C°°([0, R))and the couple (u,v) satisfies the ODE system 
(11.91) and the boundary conditions w(0) = d'(0) = 0, u(R) = ~^ s + and v(R) = —^s + . 

Conversely, if u G Hf oc ([0, R)-, rdr) fl T 2 OC ([0, R)\ y) and v G /// oc ([0, R):rdr) satisfy 
(11.9 p with the boundary condition u(R) = -^s+ and v(R) = —then the tensor Q = 
v(r)E 0 + u(r)Ei(ip) belongs to H^JBr, .5^) and is a k-radially symmetric solution of (II.4p 
and 02 .bp . 

Proof. Assume that Q G H( oc (B R , CPq) be a fc-radially symmetric solution of (II.4p and (12.bp . 
Then Q can be expressed in the form (12.51) . Standard elliptic regularity implies interior 
smoothness of any solution Q G H( oc (B R ) of (11.41) (see for instance [26]). In particular, 
Wi = Q ■ Ei are smooth on (0, R). We prove the remaining claim in several steps: 

Step 1: We prove that w\w 2 — w' 2 w 1 = 0 in (0 ,R). Fix 0 < r 0 < R. Using (12.31) and 
Corollary 12.21 for Q that is a (smooth) h-radially symmetric map in B ro , one computes that 

9 V Q = ^(S 0 Q ~ QS 0 ) in B ro , 

where Sq = 62 ® e\ — e\ < 8 > 62 - Considering now the scalar product of SoQ — QSo with both 
parts of (11.41) . we obtain 

AQ • (SoQ - QS 0 ) = 0 in B ro . 
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Integrating by parts over the ball B ro leads to 

0 = [ AQ-{S 0 Q-QS 0 )dx= [ d r Q-(S 0 Q-QS 0 )d'H 1 . 

J Brg J dB r(j 

Using the above expression of d v Q, we deduce 

2 IT 

d r Q(r 0 , ip) ■ d v Q(r 0l ip) dip = 0, for every r 0 G (0, R). 

Combining with (12. 3 1) and (12.51) . we conclude with Step 1. 

Step 2: We prove that tc 2 = 0 in (0, R). First, note that the boundary conditions on Q read 
as wo(R) = — w\(R) = ^ and w 2 (R) = 0 (which are understood as limits if R — oo). 
Next, we show that there exists 0 < R\ < R such that w 2 (r) = 0 for all r G (R\, R). Indeed, 
since w\ is continuous and w±(R) > 0, there exists an interval (Ri,R) such that w\ > 0 
on (Ri,R). The equality in Step 1 implies that ^ is constant on (R 1: R) so that w 2 = 0 
on (Ri,R). In order to prove that w 2 vanishes in the whole interval (0,i?), we write the 
Euler-Lagrange equations (1 1.4 1) within the decomposition (I2.5j) : 

+ ~fJ^ W 1 + W 2); 

2 b 2 \ 

+ T6 w T 

2 b 2 \ 


b 2 


i =0 


w’q + -^ = w 0 (—a 2 + c 2 w i ~ 

+c 2 E ^ 2 

i =0 

+c 2 E », 2 


„ , w[ _ k 

r r 


w 1 + — —w i = w i 


n w> 2 k 2 
W 2 -\ - -w 2 = W 2 


i =0 



in (0, R) where we used that Q 2 — ^-/ 3 = w ° ^ — E 0 — ^J^w 0 (wiEi + w 2 E 2 ). Then we 
apply the Cauchy-Lipschitz theory for the 2nd order ODE in w 2 (with smooth coefficients in 
(0, R)): since w 2 vanishes in (i?i, R), we deduce that w 2 = 0 is the unique solution in (0, R). 
Therefore, Q = w 0 (r)E 0 + w 1 (r)if 1 in B R and ( Wi,w 0 ) satisfies the system (11.91) . 

Step 3: We prove tCo(0) = 0 and the regularity of w 0 . Since Q is smooth in B R , we obtain 
that Wo = Q ■ E 0 is smooth in B R . In particular, Wo extends to an even (smooth) function 
on (—R,R). Therefore w 0 G C°°([0,i?)) and Wq(0) = 0. 

Step f: We prove that ici(0) = 0 and the regularity of w\. By Corollary 12.21 we know that 
Wi G // 1 ((0, R)] r dr) 0 L 2 ((0 ,R)] \ dr^. Then w\ is continuous on (0, R) and we have for 
n,r 2 G (0 ,R)\ 


w\(r 2 ) - Wi(ri)| 


f*r 2 


Wiw[ dr 


'r 1 


< 2 




1/2 
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Since the right hand side converges to zero as \r2 — rf\ —> 0, it follows that w i is continnous 
np to r = 0. Combined again with w\ G L 2 ((0, i?); - dr), we conclude that uq(0) = 0. 

For the regularity of w i, note that w\ satisfies 

w'l + -y ~ = w i(r)g(r), r G (0, R ) (2.7) 


where g is a continuous function in [0,17). Then we have (see [PS, Proposition 2.2]) that the 
function 


CM 


w±(r) 

r \k\ 


is continuously differentiable up to r = 0 with vanishing derivative CM) = 0 at the origin. 
This implies that — — i s continuous in [0,17). Returning to equation (12.711 . we deduce 
that w\ G C 2 ([0, R)). 

Conversely, assume that u G Hl oc ([0, R);rdr) D L 2 OC ([0, R); y) and v G 77/ oc ([0, /?,); rd.r) 
satisfy (jl.9[) with the boundary condition u(R) = -j= s + and v(R) = —^s + . Then Q belongs 
to Hl oc (B R , y 0 ) (thanks to (I2.4p l and satisfies (12.611 . The system (jl.9(1 implies that Q satisfies 
(II. 4p in Br\{ 0} (see Step 2 above). Since Q G H1 oc (Br, and a point has zero Newtonian 
capacity in two dimensions, Q satisfies (ll.4p in Br. This finishes the proof. □ 

Proof of Proposition \1.3[ It is a consequence of the above result. □ 


3 Study of the ODE system on finite domains 


In this section we first show the existence of a smooth solution ( u , v) of the system (1 1.9 p on 
a finite domain (0, R) with k G Z \ {0} and a 2 , b 2 , c 2 > 0 with the boundary conditions 


u(0) = 0, z/(0) = 0, ti(fi) = v(R ) - 


(3,1) 


This solution (u,v) has a sign invariance: u > 0 and v < 0 in (0,77). Second, we study the 
qualitative properties and provide appropriate upper and lower bounds on the constructed 
solution (u,v). These bounds will be extensively used in the next section when proving 
existence of the solution on the infinite domain. 


3.1 Existence of solutions with u > 0 and v < 0 

Let R G (0, oo), k ^ 0 and a 2 ,b 2 ,c 2 > 0. In order to prove existence of a solution (u,v) 
of (1 1.9 p on (0, R) satisfying (13.ip with the desired sign invariance, we will use a variational 
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approach. First, note that a solution (u,v) of the ODE system (j 1.9 [) subject to u(R) = ^ 
and v(R) = — ^ is a critical point of the reduced energy functional: H 


S(u, v) 


£r{u,v) 



{u'Y + (u ') 2 + 





rdr, 



defined on the admissible set 


S = j(u,i>) : (0, R] -> R 2 -fin', s/rv', f=, ~fiv e i 2 (0, R), u(R) = ^=,v(R) = -^= j . 

(3.3) 

If (u,v) G SF, then u is continuous on [0, R] with w(0) = 0 (see Step 3 in the proof of 
Proposition 12 . 311 . 


Proposition 3.1. For any 0 < R < oo and k ^ 0, there exists a smooth local minimizer 
(u,v) G 3F of $ such that (13.111 holds together with u(r) > 0 on (0, R] andv[r) < 0 on [0, R]. 
Moreover, u G C 2 ([0,i?]) with lirn r ^ 0 ypy > 0, v G C' 2 ([0,i?]) and 


y/3v(r) +u(r) < 0 , Vr G [0 ,R). 


(3.4) 


Proof. We divide the proof in several steps: 
Step 1: Existence of minimizers of S’ on 


3T_ := {(u,u) G 3T : v < 0}. (3.5) 

First, we know that S(u,v) is continuous and coercive in the convex closed set 3F_ endowed 
with the strong topology ^// 1 ((0, R); r dr) n L 2 ((0, R ); ~ dr) ^ x iF 1 ((0, R)-,rdr) (see Foot¬ 
note [2 ]) . Then the direct method of calculus of variations implies the existence of a global 
minimizer (u,v) of S on the subset 3F_. The couple (u,v) satishes 


„ u' k 2 u 
u + — - — = u 


-a 2 + 



? 


v" + —>v 


-a 2 - —b 2 v + c 2 (u 2 + v 2 ) 

v 6 


2„.2 


H— -j=b-u 


distributionally in (0, R) 


(3.6) 


2 If k F 0 and R < oo, we have that S{u,v) < oo if and only if v £ H 1 ((0, f£); r dr) and u € 

H 1 ((0, R); r dr) D L 2 ((0, R ); - dr). This is due to standard Sobolev embeddings and the fact that the bulk 
energy density is bounded from below (which can be seen from the inequality \v{v 2 — 3m 2 )| < (it 2 + i; 2 ) 3 / 2 
for any £l). 
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with boundary conditions 


w(0) = 0, u(R ) 



v(R) 



(3.7) 


Since u and v are continuous in (0, i?], we have by (13. 6 p that u G C 2 ((0, R]). 

Since the energy S is invariant with respect to a change of sign of u and u(R) > 0, we 
deduce that (|u|,u) is also a global minimizer of S over 27L. The strong maximum principle 
applied to the first equation in (13.61) (for (|u|,r;)) implies 


u\ > 0 in (0, R) 


since u(R) > 0. Hence 

u > 0 in (0, R). 

Also note that on the open set \v < 0), the inequality in (13.61) becomes equality and therefore, 

u,uGC' oo ({u<0}n(0,i?)). 

Step 2: We show that limsup r _> 0 u(r) < 0. Assume by contradiction that limsup r _^ 0 v(r) = 0. 
By construction, we know that v is a global minimizer of S(u, •) over the set 7_. Note now 
that in S(u,v), the contribution of v to the bulk potential is 


f^ v \r, v ) = r 


b 2 2 . —a 2 + c 2 u 2 (r) 2 b 2 3 c 2 4 

—=u 2 (r)v H- -v 2 - -=v 6 H- v 4 

VQ 2 3^6 4 


r G (0,i2),u < 0. (3.8) 


Since u(r) > 0 for all r G (0, R], we deduce the existence of 5 G (—0) such that f(r, •) 
is increasing in [5, 0] for every r G (0, R). (We highlight that 5 depends only on a 2 ,b 2 ,c 2 > 0 
and IMlicc and 5 is independent in r > 0 due to the form of the linear and quadratic terms 
in /^)(r, v).) By the above assumption, there exists an interval (Ri,R 2 ) C (0, R) such that 
v(R 2 ) — S, v > S in (i?i, R 2 ) and either R± = 0 or v(R\) = 5. Set v = v in (0, R\) U (R 2 , R) 
and v = 6 in (Ri,R 2 ). Then S(u,v) < S(u,v ) which contradicts the minimality of v. 

Step 3: We prove the following result: Let (u, v ) be a solution of (13.61) and (13.71) such that 
u > 0 and v < 0 in (0 ,R). Provided that limsup r ^. 0 v(r) < 0, then (13.41) holds true (which 
implies v < 0). Consequently, if (u,v) is a minimizer of S in then (u,v) is a local 
minimizer of S in 27. First, we define the function 


v 1 

W = - -\ - 7 = 

U y/3 


in (0, R). 


Then one computes that 


w 


/ 


v' u'v 

O 5 

u u z 


v" u"v 2 u'v' 2 {u') 2 v 

u u 2 u 2 u 3 


(3.9) 
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that leads to 


„ w 2 u If „ v 

w H-1- w = - (v + - 

r u u \ r 


v / „,u 

— ( u d 

u z V r 


(3.10) 


Using the ODE system (13.61) in (13. 10p . we obtain 


.. . 1 2w'\ . 36 2 m 

«/'+ - +— W-^ 

r «/ v/6 m V3 


v 1 

- + —p) W > 


& 2 u 

r 2 n 


in (0, /?,). 


(3.11) 


—v— 

>0 


>0 


By definition of w, since limsup r ^. 0 v(r) < 0, ii(0) = 0 and u > 0 in (0,i?), we have w < 0 
in a neighborhood of 0. By (13.71) . we also have w(R) = 0. Applying the strong maximum 
principle for (13.lOj) on (0, R ), we deduce that w < 0 on (0, R) and (13. 4 j) is now proved. 

Step 4 : We prove the regularity of u, v and lim r _> 0 > 0. By Proposition 12.31 the tensor 
Q defined by (11.81) is a smooth fc-radially symmetric solution of (II .4)1 . and so by the same 
proposition, u,v G C 2 ([0,i?]) and w(0) = n'(0) = 0. Now, let £ = -f^, then £ is continuous 
up to the origin (see the paragraph following (12.71) 1 and satisfies 

£" + (1 + 2 \k\)— -g(r)C(r) = 0 in (0, R) 
r 

for some function g G (7([0, R)). Applying [19], Lemma B.2], we see that £(0) > 0, as 
desired. □ 

Open problem 3.2. By construction, the solution (u,v) in Proposition \3. 1\ is a local mini- 
mizer of S’ over Sf. Is (u,v) a global minimizer? 

Remark 3.3. The following upper bound and uniqueness result are standard and hold for 
any a 2 , b 2 , c 2 > 0: 

1. If R G (0, oo) and (u,v) is a solution of the ODE system (11.9)1 subject to (13. ip . then 
the following upper bound holds: 

u 2j rv 2 <^-s 2 + m (0, R); (3.12) 

o 

this remains true even if b 2 = 0 (see e.g. Proposition 3 in l2Bj ). 

2. There exists R 0 > 0 (depending on a 2 , b 2 , c 2 ) such that for any R G (0, R 0 ), there exists 
a unique solution (u,v) of the ODE system (11.91) with (13.11) . This is a consequence 
of the Poincare inequality (see for instance in the related Ginzburq-Landau framework 
Thm. VIII.7, p. 98, |f;. 
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Figure 1: Schematic graphs of u and v in different regimes of a 2 , b 2 and c 2 . 


3.2 Upper and lower bounds for (w, v) 

Now we are ready to prove upper and lower bounds for any solution (u, v ) of the ODE 
system (11.91) - (13.11) with u > 0 and v < 0. These properties will be essential in proving the 
convergence of solutions on bounded domains to a solution on infinite domain. It turns out 
that these bounds strongly depend on the relation between material parameters a 2 , b 2 and 
c 2 . In fact, we will distinguish regimes leading to different behavior of v (see Figured]): 

• If b 4 > 3a 2 c 2 then v > — ^=. 

• If b 4 = 3a 2 c 2 then v = — 

• If 6 4 < 3a 2 c 2 then v < — 

The regime b 4 = 3 a 2 c 2 can be considered as a special case the other regimes. However, it 
has a distinctive feature that v = — is a local minimum of the v- relevant part of the bulk 

energy density (i.e. the function /w defined in (I3.8|) ). This allows us to establish stronger 
statement, for example the uniqueness result in Proposition 13.51 below. 

3.2.1 The regime b 4 > 3a 2 c 2 
Throughout this subsection we always assume 

b 4 > 3a 2 c 2 . (3.13) 
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Under this assumption, the following inequalities hold (see (II. 3ft ) 


where 


s_ = 


> s + > 62 
y/6~ V^C 2 ’ 

(3.14) 

b 2 — y/b A + 24a 2 c 2 

- r~9 - < o. 

4c 2 

(3.15) 


When the inequality in (13.13ft is strict, the inequalities in (13.14)1 are also strict. 

We prove the following bounds on u and v. 

Proposition 3.4. Assume b A > 3a 2 c 2 > 0, 0 < R < oo, k ^ 0 and let (u,v) G ST be any 
solution of (11.91) - f!3. Ill with u > 0 and v < 0 in (0, R). Then 


and 


^=< v <f 2 -s.m(0,R) 


uj < u < —A in (0, R ) 

V 2 


(3.16) 

(3.17) 


where ui(r ) is the unique solution of the following problem 
u'j k 2 2 b 2 V2 4c 2 


v'l + ~ ~ = u i 


—a 2 — 


Uj H——w 2 , M/(0) = 0, to(-R) = —L 

\/2 


(3.18) 


3 3 

Proof. We divide the proof in several steps: 

Step 1: We prove the upper bound v < i/1s_ in (0,R). Assume by contradiction that 


the maximum of v is attained at some point r 0 G [0 , R) where 0 > v(r 0 ) > y ^s_. By 
Proposition 13.11 we apply the maximal principle for the PDE satisfied by v G C 2 ( B /,<): 

b 2 v(r 0 ) 


Au(r 0 ) = n(r 0 ) -a 2 -pA- + c 2 y 2 ( ro ) + l v(r 0 )c 2 + ^ u 2 (r 0 ) 


<o 


V 6 


V6J 


-NT— 

>0 


>0 


which leads to a contradiction. 

Step 2: We prove a weaker lower bound 

v(r ) > 


6 2 


\/6c 2 


in (0, R). 


(3.19) 
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Assume by contradiction that the minimum of v is achieved at rq e [0, R) with v{r\) < 


-^ 3 . Then, as at Step 1, we obtain 


A v(ri) = v(ri) —a 2 — 1 + (?v 2 (r 1 ) + ( v(r 0 )c 2 + A_ ] n 2 (r 0 ) 


>o 


76 


76 


—V"* 

<0 


<0 


which leads to a contradiction. 

Step 3: We prove the optimal lower bound v(r) > —^ in (0 ,R). Using (13.4p and (I3.19P we 
obtain 


Av = v 

< v 

< V 


2 b 2 V . 2 2 
—a- ■= + c v 

76 

2 b2v . 2 2 
—a-= — c v 

7e 

_ a 2 + 2 ^ + 4cV 

76 


+ I uc 2 + u 2 


76/ 

+1 ^ + 7) 3 " 2 


in £> 


R- 


Applying the maximum principle as at Step 2 , we obtain the desired lower bound. 

Step 4 : We prove u(r) < in (0 ,R). Indeed, this upper bound follows directly from 
inequalities (13.41) and (13.161) . 

Step 5: We prove the lower bound of u. By (13.4p and (13.191) . we have 


; 2 b 2 

v ~T3 u - 2v -~7& ill(0 " R) 

Multiplying with v + < 0, we obtain: 


b z v + c 2 v 2 < - 


c 2 u 2 


3 \/3 3 


in (0, R). 


By (11.91) . the last inequality implies that u is a super-solution for (13.18p . i.e., 


// u ' k 2 / 

u -|- -u < u 


b 2 72 4c 2 2 , . 

—w + ) m (0 ,R). 


By [16] (see also n). we know that there exists a unique solution m of (13.18ft that satishes 
0 < uj < in (0, R)\ moreover, by [TjJl Proposition 2.2], we have that uj{r) = ar^ + o(l) 
as r —* 0 for some a > 0. By Step 4 in the proof of Proposition 13.11 we deduce that 
u(r) = ar^ + o(l) as r —* 0 with a > 0, we can apply the comparison principle (see [I9l 
Proposition 3.5]) to obtain that u > Uj in (0, R). □ 
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When b 4 = 3a 2 c 2 , we have the following stronger result: 

Proposition 3.5. Assume that b 4 = 3a 2 c 2 > 0, k ^ 0, R G (0, oo) and let (u,v) G PA be 
any solution of (1 1. 91) - (13. 1 1) with u > 0 and v < 0 in (0 ,R). Then v = — and u is the 
unique solution Un of the following problem: 


U "u + — - = c 2 u n {u 2 n - ^), in (0, R), 

r r z Z 

uji(O) = 0, ujj(R) = — i. 

Moreover, ujj is an increasing function with 0 < ujj < ^ on (0, R ). 


(3.20) 


Proof. Note that y|s_ = —^ Thus, by Proposition 13.41 (namely (I3.16p h v = 

— Substituting this in (II.9]) . we obtain that u satisfies the problem (13.20(1 . By [16] (see 
also m Theorem 1.3]), we know that the problem (I3.20P admits a unique solution un- 
Moreover, 0 < ujj < ^ on (0, R ). □ 

Remark 3.6. In the case R = oo, we recall that problem (I3.20p has a unique solution un in 
(0, oo) (see fT9 1 Proposition 2.5]) and the behaviour of ujj at infinity is given by: 


un(r ) 


5 + 

V2 


P 


+ o(r 2 ) , 


as r —$■ oo, 


where f3 = ^. 


3.2.2 The regime b A < 3a 2 c 2 
In this subsection we always assume 

b A < 3a 2 c 2 . (3.21) 


Under this assumption, the following inequalities hold (see (1 1.3 p 1: 


[2 s + b 2 

V 3 s - < "71 < ~ 7 ^' 


(3.22) 


We prove the following bounds on u and v. 

Proposition 3.7. Assume 0 <b A < 3a 2 c 2 , 0 < R < oo, k ^ 0 and let (u,v) G IP be any 
solution of fll.9p - fl3.ip with u > 0 and v < 0 in (0, R). Then 



< v < 


,£± 

Vq 


(3.23) 
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where uni '■ (0, R) —$■ M is the unique solution of 


n , u 'lII k2 UlII 
U IU i - 


S+ 


r2 = ncu ni (u ni - —), u HI (0) = 0, u ni (R) = 


(3.25) 


Proof. We follow the steps in the proof of Proposition 13.41 

Step 1: We prove the lower bound v > w|s_ in (0 ,R). Assume by contradiction that the 


minimum of v is achieved at some point r 0 E [0, R) with v(r 0 ) < y |s_. Then by (I3.22H . the 
PDE satished by v implies: 

= v{r 0 )[-a 2 --^=4 + cV(r 0 )] + ( v(r 0 )c 2 + ~^=)u 2 (r 0 ) 

>o vX X 


—— 

.0 


<0 


which leads to a contradiction. 

Step 2: We prove the weaker upper bound v(r) < i n (0, R). Assume by contradiction 

that the maximum of v is achieved at r\ E [0 ,R) with v(r\) > —Similarly, we obtain 


Au(ri) = v(rf)[—a 2 - ^ + c 2 v 2 [r 1 )] + (v(r 0 )c 2 + — =)w 2 (r 0 ) 


<o 


V6 


Vq j 


-V- 

>0 


-V— 

>0 


which leads to a contradiction. 

Step 3: We prove the optimal upper bound v < — in (0, R). By (13.41) . Step 2 leads to 

Av = v[—a 2 — 4L + c 2 v 2 ] + ( vc 2 + —j=)u 2 
y6 V6 

> v[—a 2 — 4= + c 2 v 2 } + (vc 2 + 4_)3 V 2 

y6 \/6 

> v[—a 2 + + 4cV]. 

VE 


As above, the maximum principle yields the desired upper bound. 























Step 4-' We prove the upper bound u < in (0, R). The inequality u < in (0, R) follows 

directly from (13.12j) and Step 3. Also, by (I3.23p . 

2 h 2 1 

—a H- -= v + c v > — -c s, 

^6 - 2 + 


and so u satisfies 


u" H —u --a > c 2 u(u 2 -si 

r r z 2 


Thus, as u < ^ and w(0) = 0, the strong maximum principle implies that u < ^ in (0, R ). 
57ep 5: IFe prove i/ie /ower bound of u. First, note that (13. 12)1 yields 

— a 2 H— -=v + c 2 {u 2 + r? 2 ) < —a 2 H— -r + /ic 2 (w 2 + v 2 ) + (1 — fi)c 2 -s 2 l _ : (3.26) 

V6 v6 3 

where 0 < p < 1 will be chosen so that the function 

£ 0 ) = ~^ v+ ^ v2 

is maximized on [<y/|s_, —^=] at the point —For that, we need to insure that ^=) > 
s_) which is equivalent to 
2 b 2 p 2 s 

—= + /rc 2 (7-s_- j=) > 0 i.e., (8 + p)b 2 — 3pVb 4 + 24a 2 c 2 > 0 

76 V 3 -v/6 


(here we used (j 1.3 [1 and (13.15)1 ). Thus, the choice /a = e (0,1) fulhlls our objective 

and we conclude that 

£(v(r)) < f (--^±=) = ~J S + + r e (°> R )- 


Combined with (13.261) . we obtain: 
2 b 2 


b 2 


a 2 + —=v + c 2 (w 2 + v 2 ) < -a 2 + /ic 2 w 2 - ^-s+ + ^-s+ + (1 - /x)c 2 ^s;7 
V6 3 6 3 


(3.27) 


Together with (11.31) . it yields 


2 h 2 s 2 

-a H- —v + c (u + v ) < uc (u -—) 

76 ; - r V 2 


and by (13.6)1 . 


w" + 


«' fc 2 




U < !IC 2 u[u 2 -7); r £ (0) R )- 


(3.28) 


As at Step 5 in Proposition [3741 we conclude to the desired lower bound using the comparison 
principle (see [H2 Proposition 3.5]). □ 
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4 Study of the ODE system on the infinite domain 


In this section we study the ODE system (1 1.9 j) on the infinite domain (0, oo) for k G Z\{0}. 
Using results of the previous section, we first prove the existence of a solution of (jl.9p subject 
to (ll.lOp . As consequence, we prove existence of fc-radially symmetric solutions of (II.4^ on 
the whole M 2 stated in Theorem 11.41 Next, we prove finer asymptotic behavior at infinity of 
any solution of (II ,9jl subject to (11.101) . 


We start by proving the following existence result on (0, oo). 

Proposition 4.1. Let a 2 , b 2 , c 2 > 0 be fixed constants and k G 7L \ {0}. Then there exists a 
smooth solutior^_ \ (u,v) of (11,91) defined on (0, oo) with boundary conditions (11.101) . Moreover, 
we have 0 < u < v < 0 in (0, oo) and (u,v) is locally minimizing in the following sense: 


<£r{u, v ) < £r(u + £, v + 7]) for any (£, rj) G C£°(0, R) with sup \rj\ < min 

(0 ,R) 

for any R > 0, where Sr is given by (13.21) . 

Proof. For every n G N*, let (u n ,v n ) be the solution of (jl.9j) on the interval (0, n) subject to 
(13. ip constructed in Proposition 13.II We extend u n and v n to the infinite domain (0, oo) by 
setting the functions (u n ,v n ) : (0, oo) —* 1R 2 as follows: 


(— 




Un, 



r ) r G (0, n ) 
r > n 


v n {r) r G (0, n) 


Since {(u n , v n )} n >i are uniformly bounded in L°°(0,oo), we have by standard regularity 
arguments that for any given compact interval J C (0, oo) and for large enough n > nj, the 
couples {(«„, v n )} n >nj are uniformly bounded in C 3 4 (J). Using the Arzela-Ascoli theorem, 
we deduce that u n —>■ u and v n —> v in C) 2 c (0, oo) (up to a subsequence). Thus, (u,v) : 
(0, oo) —> M 2 satisfy (1 1.9 p on (0, oo), too. Bv Propositions 13.5113.4113.71 Step 3 in the proof 
of Proposition 13.11 and Remark 13.31 we hav([j 


min 






< v < max 


u + V3v < 0, 


2 I 2 

U + V 


< 








3 Here, u and v are C 2 up to r = 0. 

4 For the strict upper bound of u, see Steps 4 in the proofs of Propositions 13.41 and 13. 71 
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Also, u > 0 in (0, oo). 

We next show that u G -£// oc ([0, oo); r dr) 0 L 2 ([0, oo); ^) and v G Hl oc ([0, oo); r dr). 
Thanks to the (uniform) bound of u n and v n , it suffices to show that <o m (u n) v n ) is uniformly 
bounded for n > m > 0. Indeed, if we (u m ^ n (r), Vn hn {r)) be the extension of ( u m ,v m ) which 
equals to ( u n ,v n ) in the interval (m + l,n) and is linear in [ m,m + 1], then 


0 A ^n(^n, ^n) X m,n ) 

( ?m(V'ni Un) ^rn( t ^rn ) 


+ 



fc 2 

(|^ n | A |u n | "h y |Wn| ) "P fbulk(u n , U n ) vdv 
£ 2 

(I ^m,n I + km,nl +“d ^m,n | ) + fbulk(u ra,n? ^ m,n ) 


rdr , 


where (by a slight abuse of notation) 


fbuik(x,y) = fbuik{xE 1 + r/£ 0 ) = -y(a ; 2 + 2 / 2 ) + y (^ 2 + Z / 2 ) 2 - y^y 2 /(d / 2 “ 3a;2 )- 

As ( u n ,v n ) are uniformly bounded in (m,m + l), (u m:n (r), v m)7l (r)) and its derivative are also 
uniformly bounded in ( m,m + 1). It thus follows that 0 > £m{u n ,v n ) — 4(«m,t'ra) — C for 
some constant C independent of n. This proves that ^ m {u ni v n ) is uniformly bounded. 

The locally minimizing property of ( u , v) follows from the bounds for v n and the minimiz¬ 
ing property of (u n ,v n ). It remains to show that (u,v) takes on the desired value at infinity 
(the boundary condition at the origin and the smoothness of u and v are a consequence of 
Proposition 12.31) . 

Case 1: b 4 = 3a 2 c 2 > 0. By Proposition 13.51 we have that v = — Moreover, since u n is 
the unique solution of (I3.20p in (0,n), we know by [US Proposition 2.4] that u n converges in 

Cf oc to the unique solution u = un of (I3.20p in (0, oo), and so w(oo) = 

Case 2: b 4 > 3a 2 c 2 > 0. By Proposition 13.41 the same argument as above implies that 

— < v < \J\ S - < 0 and 0 < Ui < u < in (0, oo) where Uj is the unique solution 

of (13.181) in (0,oo), in particular, zt(oo) = For n(oo), note on one hand that v > 

— ^ in (0, oo) so that liminf r ^. 00 n > — y. On the other hand, by (13.4p . we have that 

limsup^^u < — ^lim^oott = — We thus have u(oo) = — ^=. 

Case 3: 0 < b 4 < 3a 2 c 2 . Arguing as in the previous case (but using Proposition 13.71 instead 
of Proposition I3f4l) . we get w(oo) = ^t. Next, since v < —^ in (0, oo), lirnsup r ^ 00 v < — %. 
On the other hand, by (13.121) . 


lim inf v > — 

r—>■ oo 



w(oo) 2 


g+ 

VE' 
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We again obtain v(oo) = — ^ as desired. □ 

We now prove the existence of fc-radially symmetric solutions of HI.4)1 subject to ('11.511 : 

Proof of Theorem l.f. The assertion is a consequence of Propositions [T73l and |4~T1 □ 

In the proof of the instability result, we need some detailed behavior at oo of any solution 
(u,v) of the system (11.911 subject to (11.1011 : 

Lemma 4.2. Let u mid, v be any solution of (II.9H defined on (0, oo) subject to (II . 10)1 . Then 
( u,v ) has the following behavior as r oo: 


u = 


s + V 2 k ' 2 2 b 2 + c 2 s + 


72 
s+ 




2 b 2 (— b 2 + 4c 2 s + ) 
a/6 k 2 —b 2 + c 2 s + 


r 2 + 0 (r 4 ), 


r 2 + 0 (r ). 


a/6 2 b 2 (-b 2 + 4c 2 s+) 

The proof of this result uses the following lemma: 

Lemma 4.3. Let Br c R n with 0 < R < oo. Assume for some constant C > 1 that 


(4.1) 

(4.2) 


— < h(x) < C inM. 71 \ Br. 
o 

If u G C 2 (M n \ Br) satisfies 

—A u + h(x ) u = 0 (\x\~ a ) 

for some a > 0 and if u{x ) —> 0 as |x| -A oo, then u = 0(\x\~ a ), where the big“ O” notation 
is meant for large |x|. 

Proof. Let L = — A + h(x). We have 

L(\x\~ a ) =a(a-n + 2)\x\~ a ~ 2 + h(x) |x|- Q . 


Hence, by our assumption on h(x), 

\x\~ a < L(\x\~ a ) < C \x\~ a in W 1 \ B 2 r. 

o 

It thus follows that, there is some large radius R' > 2 R and some C\ > 0 such that 
L(u - C^x]-*) < 0 < L{u + C\ |x|- Q ) in M n \ B R ,. 

Replacing Ci by a larger constant if necessary, we can also assume that 

u — Ci\x\~ a <0 < u + Ci\x\~ a on DBr'. 

The assertion follows from the maximum principle. □ 
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Proof of Lemma f.2 


Let u — u — u{ oo) and v — v — v(oo). We have 


u" + ^u' = (ci + 0(\u\ + \v\+r 2 ) j u+ (c 2 + 0(\u\ + \v\v + (4.3) 

v" + -v' = [c 2 + 0(\u\ + |0|)^u + (c 3 + 0(\u\ + |u|))t). (4.4) 

where c\ = c 2 s+, c 2 = — ^(c 2 s+ — 6 2 s + ), and c 3 = |(26 2 s + + c 2 s 2 ). 

Introducing X — u + \/3h and Y = \/3w — L, we obtain 

X" + l -X' = (b 2 S+ + 0(1*1 + |u| + r“ 2 ))x + 

Y" + -Y' = - U c 2 s 2 - b 2 s+ + 0(|-u| + |*| + r" 2 )) Y + g+ . 

r 3 V / V2r 2 

Since both 6 2 s + and 4c 2 — 6 2 s + are positive and since 0(|*| + |h|) = o(l) as r —>■ oo, 

Lemma [4.31 implies that |X| < C r -2 and |Y| < Cr~ 2 . It follows that the above equations 
of X and Y can be rewritten as 


1 p Q 

x" + -X' = b 2 8+ X + + 0(r~ 4 ), 

r V2r 2 

Y" + -Y' = i(4c 2 4 - 6 2 S+ ) y + + 0(0. 

r 3 \/2r 2 


Thus, the functions X = X + r 2 and Y = Y + 2 satisfy 

X" + -X' = 6 2 s + X + 0(r -4 ), 
r 

Y" + ly' = -(4c 2 s 2 - b 2 s+) Y + 0(r~ A ). 
r 3 

Again, Lemma [4.31 implies that |X| + |Y| < Cr~ A . Returning to the variables u and v, we 
obtain the desired asymptotic expansion. □ 


5 Instability of /c-radially symmetric solutions 

In this section we prove the instability of radially fc-symmetric solutions of (11,4[) on the whole 
M 2 for \k\ > 1. Note that for any Q <E il/ oc (K 2 , y 0 ) satisfying (11.5ft . one has &(Q) = oo. We 
thus adopt a second variation at Q in a local sense as defined in (11.111) . 
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Proof of Theorem \1.7[ We follow the ideas from |18| 20]. For \k\ > 1, let Q be a fc-radially 
symmetric solution of (jl,4[) on M 2 subjected to (11.51) . Then Q has the form (11. 8 ft with (u,v) 
satisfying ( 11 . 9|) - ( 11 . 1011 . 

Let £ > 0 be a small parameter. Since w(oo) = there exists R > 0 such that 

(1 — e)u(oo) < u < (1 + £)w(oo) in (P, oo). (5.1) 


We take 


P = w(r)h((p)—=(n ® e 3 + e 3 <g) n) 
v 2 


where n is as defined in (II. 61) . w G Cf°(R, oo) and 


Then P e C?(R 2 ,y 0 ). 
We have 



if k is odd, 
if k is even. 


dP 1 dP 1 k 

—— = w'h—= (n <g) e 3 + e 3 <g) n ), tt - = wh! —= (n <g) e 3 + e 3 <g> n)+wh — -= (m <g) e 3 + e 3 <g) m) 
or V 2 cty? \/2 2 V 2 


so 


|VP | 2 _ |<9 r P | 2 |a v P | 2 _ (w/h ) 2 w 2 (//) 2 w 2 A: 2 h 2 

2 ~~ 2 + 2r 2 ~~ 2 + + 2r 2 + 8 r 2 

We also have: 


-2tt 


-2?r 


h 2 (ip)dip = 7 T, / (ti((p)) 2 dip = c k = -(1 + (-l) fc+1 ), 


'0 J 0 

}2 1 .„2l2 1,012 „,2 i „,2 


|P 2 | = w 2 h 2 , \Q\ 2 = u 2 + v 2 , tr(PQ) = 0 


o w 2 h 2 . . . 9 . w 2 h 2 . u v . 

P = — (n ® n + e 3 ® e 3 ), tr(P Q) = —(— + —), 


hence (II . 11 jl becomes: 

-JS?[Q](^) = f 

71 Jo 


k T Cfc 9 / 2, 2 2 2 / 2 i 2 \ \ 2 


k '| 2 + 


4r2 —+ ^ 2 r> + c\u 2 + n 2 ) ) w 2 


b 2 


^+ Vsv'j tc 2 |r dr, 
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where 


Cfc = i(l + (-l) t+1 ). 

We now use the Hardy decomposition trick as in [IE, 120] by setting w = uf with £ £ 
Cf°(R,oo). Then: 


7T 


•^[Q\(P) — / \ l^l +1^1 + 2uu'ff' + 


, k 2 + c k 


u 


Ak 2 u 2 


4r2 u*e+(u"+--- 4r2 


K 2 


^-= (u + VSv'j w 2 £ 2 jr dr , 


V2 


le 


/12 


3 k 2 - c fc 2 


b 2 


4 r 2 + 


where for the first equality we used the equation (11.91) for u and for the second equality we 
integrated by parts the term f Q u"uf 2 rdr. From Lemma f4~2l we know 


u 


V3v = 


k 2 


V2 b 


:2 r 2 


+ 0(r ) as r —> oo. 


Therefore, by replacing R by a larger constant if necessary, we can assume that 

b 2 / rr \ —life 2 + c k 


V2 


(u + V3v j > 


20r 2 


in (R, oo). 


Hence, for any £ £ Cf°(R, oo), we deduce by (15. lj) and k 2 — c k > 3: 


~^mp) < 


71 


'R 


le 


l\2 


k 2 - c k 
5 r 2 


£ 2 > u 2 r dr < (1 + £) 2 w(oo) 2 


' R 


I £\ 2 -^e\rdr. 


It is not difficult to find a test function £ 0 £ Cf°(R, oo) such that El 

1 


IR 


l£ol 2 


2r 2 


£o } r dr < 0. 


The result follows immediately. 


□ 


Remark 5.1. Theorem 23 and its proof provide an insight into the stability of the k-radially 
solution on finite domains B R ( 0) for R small, respectively R large. 


For example, take £o to be a smoothing of sin(-^ 1 ) l( e xp( 2 n 7 r).exp( 2 (n+i) 7 r)) f° r some n sufficiently large. 
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• Case 1: if R is small (ie, R < Ro(a 2 , b 2 , c 2 )) then one can use the Poincare inequality 
to show that the solution of the PDE system (11.41) with boundary conditions (I2.6[) is 
unique (see for instance in the related Ginzburg-Landau framework Thm. VIIL7, p. 

This unique solution must necessarily be the global minimizer of ^ onfl = Br 
subject to (12. 6 j) . Therefore, it coincides with the k-radially solution of Proposition 2.3 
(write HO) and it is stable as global minimizer of TP. 


Case 2: if R is large (ie, R > Ri(a 2 ,b 2 ,c 2 )) then for \k\ > 1 the k-radially solutions 
obtained in Proposition ^. 1\ are expected to be unstable, since the solutions ( ur,vr ) of 
the ODE system on the finite domain [0, R] will suitably approximate the solution in 
the whole space obtained in Proposition^ 7j 


A Remarks on the case b 2 = 0 

In this appendix, we collect some known results regarding the case b 2 = 0 and k e Z \ {0}. 
In |3], it was shown that on a Unite disk the system (jl.9[) and (13.ip has a unique solution 
(u,v) with the sign invariance u > 0 and v < 0. Furthermore, Q = u Ei + v E 0 is the 
unique global minimizer of the full Landau-de Gennes energy Jp subjected to the boundary 
condition (12.61) . 

For infinite domain, the situation is different. We have: 

Theorem A.l. Assume that b 2 = 0, a 2 > 0, c 2 > 0 and k ^ 0. There is no solution of the 
boundary value problem (II .9[) - (jl. 10[) which satisfies u > 0 and v < 0 in (0, oo). 

Proof. Indeed, assume by contradiction that there exists a solution ( u,v ) of (11.91) on (0, oo) 
subject to (11.101) with v < 0 in (0, oo). By (13. 12ft and (11.31) . we have that —a 2 +c 2 (u 2 +v 2 ) < 0 
in (0,oo). Hence, the equation (II.9ft for v < 0 implies (tv')' > 0 for every r > 0. Since 
i/(0) = 0, we deduce r i —> rv'{r ) is a nonnegative and nondecreasing function. It follows 
that, for any s > r > 0, v(s ) > v(r) + rv'(r) In -. Fixing r and taking a limit s —> oo it is 
clear that starting from some point So function r;(s) becomes positive. Since v is negative in 
(0, oo), this implies that v' = 0 in (0, oo). By (11.101) . we thus have v = — fA, in (0, oo). Using 
the second equation in (II.9j) . we obtain — a 2 + c 2 {u 2 + v 2 ) = 0, and so, by (II. 101) . u = % is 
constant. This contradicts the hrst equation in (11.91) . □ 
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